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USE OF SINGLE NUCLEOTIDE POLYMORPHISM IN THE CODING REGION OF 
THE PORCINE LEPTIN RECEPTOR GENE TO ENHANCE PORK PRODUCTION 

[0001] This application claims the benefit of United States provisional application serial number 
60/553,582, filed March 16, 2004, and United States provisional application serial number 
60/493,158, filed August 7, 2003 

FIELD OF THE INVENTION 

[0002] The invention relates to methods for improving swine genetics and pork production and 
to compositions and kits useful to carry out such methods and to herds produced by said 
methods. Another aspect of the invention relates to the identification and use of a single 
nucleotide polymorphism in the porcine leptin receptor (LEPR) gene. The invention is also 
drawn to the use of probes to detect the LEPR gene polymorphism in order to identify those 
animals useful for as breeding stock for improved pork production. 

BACKGROUND OF THE INVENTION 
[0003] The pork industry is experiencing phenomenal growth as it continues to meet worldwide 
consumer demand for what has become the meat product with the highest consumption. One key 
to mamtaining industry growth and cost effective production is the continued implementation of 
high-quality standards into every level of the business. 

[0004] In the United States, pork production is a vital part of the economy. Nearly 19 billion 
pounds were processed from about 97 million hogs in 2001. The economic impact of the 
industry on rural America is immense. Annual farm sales typically exceed $1 1 billion, while the 
retail value of pork sold to consumers reaches $38 billion each year. 

[0005] Pork also provides employment well beyond the farm. The U.S. pork industry is 
responsible for over $72 billion in total domestic economic activity. In addition, the pork 
industry supports over 800,000 jobs and adds over $27 billion of value to basic production inputs 
such as corn and soybeans. 

[0006] There are approximately 85,760 pork operations today compared to nearly three million 
in the 1950s. Farms have grown in size; nearly 80 percent of the hogs are grown on farms that 
produce 5,000 or more hogs per year. 

[0007] Major technological advancements have allowed for production to grow dramatically 
over the years. A number of innovations, including the use of genetic capabilities for higher 
reproductive efficiencies and enhanced lean muscle growth, capturing economies of size, and 
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developing animal management methods that have controlled diseases, have led to improved 
productive efficiency. In addition, U.S. pork producers are increasingly using state-of-the-art 
innovations designed to provide an environmentally efficient operation that ensures safe, high 
quality food for consumers. 

[0008] There is also an increasing consumer demand for meat products with specific qualities, 
(e.g., low fat content). This demand is fueled by accumulating evidence in the scientific 
literature that a high consumption of animal fat, especially fat with a high proportion of saturated 
fatty acids, represents a significant health hazard, including risk for cardiovascular disease. 
Another health concerns associated with high fat meats is their high cholesterol content. 
[0009] Faced with larger average farm size and consumers who seek a healthier meat product at 
a minimum cost, pork producers are continually pressed to reduce the cost of production and 
offer healthier products to stay competitive. 

[0010] One tool used in pork production is the use of genetic differences that exist among 
individual meat producing animals as well as among pig breeds. These differences can be 
exploited by breeding techniques to achieve animals with these desirable characteristics. For 
example, Chinese breeds are known for reaching puberty at an early age and for their large litter 
size. In contrast, European and American breeds are known for their greater growth rates and 
leanness. 

[0011] The occurrence of desirable traits (e.g. growth rate or muscle mass) in an animal and/or 
herd may be optimized by identifying those genes or genetic loci associated with variation in a 
particular trait of interest and increasing the incidence of the desirable allele of that gene or locus 
within a given pig population. This is necessary because the heritability for desired traits may be 
quite low. For example, heritability for litter size is around 10%-15%. Standard breeding 
methods that select individuals based upon phenotypic variations do not take into account 
genetic variability or complex gene interactions which may exist. Consequently, an improved 
approach that incorporates analysis of variation in an animal's DNA is desirable. Such a method 
provides a means for genetically evaluating animals to enable breeders to more accurately select 
those animals that not only phenotypically express desirable traits but also have the underlying 
favorable genetics. In theory, this can be accomplished by marker assisted selection. 
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[0012] RFLP analysis has been used by several groups to analyze pig DNA. Jung et al. 5 Theor. 
Appl. Genet., 77:271-274 (1989), incorporated herein by reference, discloses the use of RFLP 
(restriction fragment length polymorphisms) techniques to show genetic variability between two 
pig breeds. Polymorphisms were demonstrated for swine leukocyte antigen (SLA) Class I genes 
in these breeds. Hoganson et aL, Abstract for Annual Meeting of Midwestern Section of the 
American Society of Animal Science, Mar. 26-28, 1990, incorporated herein by reference, 
reports on the polymorphism of swine major histocompatibility complex (MHC) genes for 
Chinese pigs; these were also demonstrated by RFLP analysis. Jung et aL Animal Genetics, 
26:79-91 (1989), incorporated herein by reference, reports on RFLP analysis of SLA Class I 
genes in certain boars. The authors state that the results suggest that there may be an association 
between swine SLA/MHC Class I genes and certain production and performance traits. They 
further state that the use of SLA Class I restriction fragments, as genetic markers, may have 
potential in the future for improving pig growth performance. 

[0013] In order to exploit the advantages of a specific favorable genetic allele one must first 
identify at least one genetic marker for each desired trait. The marker(s) may be linked to a 
single gene or to a number of genes, providing additive effects. DNA markers have several 
advantages; segregation is easy to measure and is unambiguous. Moreover, DNA markers are 
co-dominant, allowing all genotypic classes to distinctly identified. Also, DNA marker 
information can be assessed at an early age (prior to expression of the phenotype of interest) and 
markers for sex-linked and sex-influenced traits can be measured in both sexes. 
[0014] The use of genetic differences in receptor genes has become a valuable marker system for 
selection. For example U.S. Pat. Nos. 5,550,024 and 5,374,526 to Rothschild et aL (each of 
which is incorporated herein by reference) disclose a polymorphism in the pig estrogen receptor 
gene that is associated with larger litter size, the disclosure of which is incorporated herein by 
reference. Another example is provided by U.S. Pat. No. 5,935,784, filed August 10, 1999, 
which discloses polymorphic markers in the pig prolactin receptor gene that are associated with 
larger litter size and overall reproductive efficiency. 

[0015] The leptin receptor (LEPR) gene encodes the leptin receptor protein, which is a cytokine 
receptor that specifically recognizes the ligand "leptin." Upon binding its ligand the leptin 
receptor initiates a cellular signal transduction cascade that ultimately produces major 



WO 2005/017204 



PCT/US2004/023050 



-4- 

physiological results, most significantly suppression of appetite. Expression variation in LEPR 
has been found in different nutritional states (Dyer et al). Reviews of known functions of leptin 
and the leptin receptor are provided in Barb et al and Tartaglia. 

[0016] The porcine LEPR gene has been localized to chromosome 6, at approximately 122 
centiMorgans (cM). Moreover, a number of DNA sequences (genomic and cDNA) for the 
porcine LEPR gene are available from the Genbank public DNA database, including: accession 
numbers: AF092422 (Ruiz-Cortez et al\ AF167719 (Hu etal), AF1S4173, AF1 84172 and 
AH009271 (Lacroix et al), AJ223163 and AJ223162 (Stratil et al), U72070 (Ernst et al), 
AF036908 (Matteri, R.L.), and U67739 (Matteri, R.L. and Carroll, J.A.), each of which are 
herein incorporated by reference. 

[0017] The murine autosomal recessive mutations obese (OB), diabetes (DB) and fatty (FA) 
were first reported in the 1960s. The phenotypes of animals homozygous for these mutations 
include severe, early-onset obesity, insulin resistance and susceptibility to diabetes. The OB 
gene has recently been cloned in human and mouse and its protein product identified as leptin. 
Subsequent research led to the identification of a receptor for leptin in mice (OBR). The gene 
for OBR was shown to map to within a 5.1 cM interval of mouse chromosome 4 that also 
contains the db locus. This report was followed by two studies providing evidence that db is the 
gene encoding OBR. A recent report by Chua and associates has confirmed that db,fa and obr 
are all mutations of the same gene. The mouse leptin receptor gene has now been assigned the 
symbol, Lepr, which replaces the previously used symbols OB-R and obr. Mapping of human 
leptin receptor gene (LEPR) has also recently been reported. 

[0018] The leptin receptor in mice (and humans) is a class-I transmembrane cytokine protein 
existing. in two forms (i.e., forms having either a short or a long cytoplasmic domain). Only the 
long form is believed to be capable of signal transduction. In mice, the LEPR gene product is 
believed to bind leptin (the 146 amino acid protein secreted into the blood by fats cells) in a 1:1 
ratio (Devos et al, Dyer et al, Tartaglia). Administration of leptin to ob/ob mice, which are 
deficient in the production of leptin, causes a reduction in food intake and weight loss (Devos et 
al). In ewes the LEPR is expressed in the anterior pituitary and adipose tissues. Moreover, it is 
differentially expressed in well-fed versus feed-restricted ewes (Dyer et al). 
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[0019] It has been hypothesized that various polymorphisms in pLEPR may affect commercially 
significant traits. For example, U.S. Pat. No. 6,458,531 (Rothschild et al\ Strait et al and 
Vincent et al (which are herein incorporated by reference) describe genetic markers, based upon 
polymorphisms in and around the pLEPR gene. These polymorphisms are described as relating 
to leanness in pigs. The Rothschild et al. £ 531 patent suggests that use of the pLEPR markers 
described therein would permit genetic typing of pigs for their pLEPR allelic variants and for 
determination of the relationship of specific RFLPs to leanness. Thus, it is suggested that the 
described markers may be used as a selection tool in breeding programs to develop lines and 
breeds that produce litters containing offspring with less fat content. 

[0020] However, none of the pLEPR polymorphisms described thus far are believed to cause any 
variance in the protein encoded by the pLEPR gene. Moreover, no determination of their nature 
(other than the fact that they are restriction fragment length polymorphisms) has been reported. 
[0021] Study of the mouse LEPR indicates that the leptin binding domain resides in amino acid 
residues 323-640. Furthermore, co-expression of the active form of the receptor with an inactive 
mutant indicates that in its functional form the receptor may exist as a multimeric complex in the 
absence of leptin (Ming et al). 

[0022] Ovilo et al have investigated the LEPR gene as possibly affecting carcass composition in 
pigs. When testing the RELP previously published by Stratil et al they confirmed an association 
between that polymorphism and fatness, but concluded that the RFLP was merely in some level 
of linkage disequilibrium with the causal mutation. The authors attempted to test the strength of 
association between carcass composition traits and the two RFLPs described in the Rothschild et 
al. '531 patent, but could not find an association because the polymorphisms did not occur 
frequently enough in the population tested. 

[0023] In view of the discussion above it is likely that both leptin and the leptin receptor product 
play some part in the determination of body composition, fatness, muscle leanness, and feed 
intake in swine. Therefore, there exists a need to identify genetic markers that are linked to the 
expression of desirable commercial traits in pigs. There is also a need for methods of identifying 
the presence of absence of these markers in individual animals and of using these markers as part 
of a pig management or pig breeding program for the improvement of pork production. The 
invention described herein satisfies these needs. 
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SUMMARY OF THE INVENTION 

[0024] To meet the needs described above, the present invention provides a method for 
screening pigs to determine those that will be likely to produce offspring with desirable genetic 
traits. The method comprises: 1) obtaining a sample of genomic DNA from a pig; and 2) 
analyzing the genomic DNA obtained in 1) to determine which pLEPR alleles(s) is/are present. 
The information collected by this method may then be used in preparing a breeding plan for 
increasing the frequency of the desired allele. 

[0025] The instant invention is further drawn to methods that comprise determining which 
variant of the pLEPR polymorphism is extant in an animal, or a plurality of animals, and then 
using this determination to formulate a breeding plan to increase the frequency of the desired 
allele and/or improve the quality of offspring produced. 

[0026] According to one embodiment of the present invention, one useful allelic polymorphism 
comprises a "C/T" variation in the fourth exon of the leptin receptor gene. This variation results 
in the pLEPR protein produced from these variants having either a methionine or a threonine as 
amino acid number 69 of the prepro pLEPR protein. 

[0027] Various embodiments of the invention provide methods for detecting which allelic 
variant is present in a particular animal. These methods include, but are not limited to DNA 
sequencing, restriction fragment length polymorphism (RFLP) analysis, heteroduplex analysis, 
single strand conformational polymorphism (SSCP), denaturing gradient gel electrophoresis 
(DGGE), polymerase chain reaction (PCR), real time PCR analysis (TAQMAN®), temperature 
gradient gel electrophoresis (TGGE), primer extension, oligo-specific hybridization and 
INVADER® genetic analysis assays. 

[0028] The INVADER® platform is based on a "perfect match" enzyme-substrate reaction. The 
INVADER® reaction uses proprietary CLEAVASE® enzymes, which recognize and cut only 
the specific structure formed during the Invader process. Instead of relying on target 
amplification, as in traditional methods, the INVADER® reaction generates its own signal 
amplification. In the INVADER® process, two short DNA probes hybridize to the target in the 
presence of the variation of interest to form the structure recognized by the CLEAVASE® 
enzyme. The enzyme then cuts one of the probes to release a short DNA flap. Each target can 
induce the release of several thousand flaps per hour. Each released flap can act as an Invader 
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oligonucleotide on a FRET (fluorescence resonance energy transfer) cassette to create another 
structure recognized by the CLEAVASE® enzyme. If recognition occurs, the CLEAVASE® 
enzyme cuts the labeled probe, which emits a detectable fluorescent signaL Each flap generates 
thousands of signals per hour, yielding millions of detectable signals per target. The 
INVADER® reaction results are easily read on most existing fluorescence detection systems. If 
the variation in question is not present, then there is no overlap with the probe, the INVADER® 
oligo, and the target DNA. Hence, there is no recognition by the CLEAVASE® enzyme, and no 
flap is released. In the absence of the cleaved flaps, no invasive structure is formed, which 
means that no fluorescent signal is released from the FRET cassette. (INVADER® and 
CLEAVASE® are registered trademarks of Third Wave Technologies Inc., Madison, 
Wisconsin). 

[0029] The current invention also provides for kits comprising the components necessary to 
carry out methods for identifying polymorphisms in the LEPR gene. These kits comprise the 
components necessary to carry out any type of analysis suitable to detect the polymorphisms 
described herein. For example, analytical methods contemplated as being useful for the instant 
invention include, but are not limited to, the following: DNA sequencing, restriction fragment 
length polymorphism (RFLP) analysis, heteroduplex analysis, single strand conformational 
polymorphism (SSCP), denaturing gradient gel electrophoresis (DGGE), polymerase chain 
reaction (PCR), real time PCR analysis (TAQMAN®), temperature gradient gel electrophoresis 
(TGGE), primer extension, oligo-specific hybridization, and INVADER® assays. Any other 
suitable means for analyzing the structure of nucleic acids is also within the scope of the present 
invention. 

[0030] Also contemplated as part of the instant invention are methods and kits for detecting the 
allelic variation at the level of protein. For example, lots comprising the components for 
immunological assays, such as the necessary antibodies, buffers, and labeling compounds, fall 
within the scope of the present invention. 

[0031] Another embodiment of the instant invention provides for the necessary novel reagents 
for the kits provided above. Various aspects of this embodiment of the invention provide for 
oligonucleotide primers suitable for use as a DNA and/or RNA probes or as primers for DNA 
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and/or RNA synthesis. Additionally, this aspect of the invention provides for antibodies useful 
for detecting proteins produced from the allelic variant provided herein. 

[0032] One embodiment provides for a method for producing pigs and the pigs produced by that 
method (considered both as individuals and as a herd). Generally, the method comprises 
analyzing either one or a plurality of pigs and determining which form(s) of the pLEPR 
polymorphism each animal possesses. Next, this information is used as part of a breeding plan to 
produce one or more pigs having the desired qualitative and/or quantitative traits. 
[0033] According to one aspect of this embodiment of the invention the method for producing 
pigs is employed to provide pigs having more desirable characteristics with respect to economic 
traits selected from, but not limited to, one or more of the following: the average feed intake 
and/or average daily weight gain, backfat, muscle mass, water holding capacity, meat color, 
intramuscular fat, meat tenderness, and/or cooking loss. 

[0034] Another embodiment of the instant invention provides a method for increasing meat 
production in a herd by a method comprising modifying the herd genetics (e.g. the frequency of a 
particular LEPR gene allele) as provided herein. One aspect of this embodiment of the invention 
provides for a method wherein the EBV of the herd is improved over time with respect to the 
trait of meat production by, for example, increasing the frequency of the LEPR gene allele which 
has been shown to be linked to increased meat production. 

DESCRIPTION OF THE FIGURE 

[0035] The following figures forms part of the present specification and is included to further 
demonstrate certain aspects of the present invention. The invention may be better understood by 
reference to this figure in combination with the detailed description of specific embodiments 
presented herein. 



Figure 


Description 


1 


DNA sequence (SEQ ID NO: 10) and Amino Acid sequence (SEQ ID NO: 11) 
of the portion of the pLEPR gene which contains the M69T and S73I 
polymorphisms. Primer sequences are underlined, the single nucleotide 
polymorphisms and accompanying amino acid changes are shown in bold. 
Nucleotide sequence without accompanying amino acid sequence is intronic. 
The forward primer starts at position 311 of Genbank accession AF1 84172, " 
Sus scrofa leptin receptor (LEPR) gene, exon 4 and partial cds". The M69T 
polymorphism is at nucleotide position 609 of AF1 84172. 
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DETAILED DESCRIPTION OF THE INVENTION 

[0036] The following definitions are provided in order to aid those skilled in the art in 
understanding the detailed description of the present invention. 

[0037] As used herein the term "average feed intake" refers to the average amount of food 
consumed by an animal or animals during a defined "period". The definition may be used with 
reference to an individual animal or, alternatively, to a group of animals such as a litter, group of 
litters, or an entire herd. The "period" may include definite time periods such as intake per day, 
week, or month. Alternatively, it may be used, for example, to refer to average feed intake 
consumed by a group of animals by the time they reach a specified stage (e.g., "weanling", 
"grower", or "finisher"). 

[0038] As used herein the term "backfat" denotes a measurement of the thickness of the fat, as 
measured on the carcass (or by ultrasound prior to slaughter), at a specified point on the animal's 
back. 

[0039] As used herein the term "BLUP" (which is an acronym for best linear unbiased 
prediction) refers to any of the various commercially available computer programs that are used 
for genetic evaluation of an animal and/or herd. Typical input parameters and data for BLUP 
programs include genetic parameter estimates, phenotypes and pedigrees. 

[0040] As used herein the term "breeding plan" preferably refers to a program for improving 
herd genetics, including the average estimated breeding value (EBV) for the herd, using the 
methods provided herein. The "breeding plan" may employ the use of statistical models and/or 
computer programs, such as BLUP, to formulate the most effective means to achieve the desired 
genetic improvement and/or allelic frequency in the herd. 

[0041] As used herein the term "cooking loss" preferably refers to the difference in weight, due 
to water loss, between a piece of meat before and after cooking.. 

[0042] As used herein the term "economic trait locus" (ETL) preferably refers to a location on a 
chromosome that is linked to "quantitative trait" providing economic value. 
[0043] As used herein the terms "efficient growth traits" and/or "performance traits" preferably 
refers to a group of traits that are related to growth rate and/or body composition of the animal. 
Examples of such traits include but are not limited to average daily gain, average daily feed 
intake, feed efficiency, back fat thickness, loin muscle area, and lean percentage. 
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[0044] As used herein the term "estimated breeding value" (EBV) preferably refers to a specific 
numeric value for an animal that predicts its "breeding value". EBV is often calculated using 
commercially available analysis programs (the output from BLUP is an example of an EBV). 
[0045] As used herein the term "fixing a genotype" preferably means producing a population of 
pigs that are all homozygous for the same allele of a particular marker in a specific gene or at a 
specific locus. 

[0046] As used herein the term "gene" refers to a sequence of DNA responsible for encoding the 
instructions for making a specific protein within a cell (including when, where, and in what 
abundance the protein is expressed). 

[0047] As used herein the term "intramuscular fat" refers to a measure of the fat content of a 
specific cut of meat (e.g. loin or ham) that is determined by chemical analysis. 
[0048] As used herein the term "linkage disequilibrium" refers to: a non-random association of 
alleles at two or more loci. A quantitative measure of linkage disequilibrium is correlated to the 
probability of two alleles (at separate loci) being inherited together. 

[0049] As used herein the term "locus" refers to a specific location on a chromosome (e.g. where 
a gene or marker is located). "Loci" is the plural of locus. 

[0050] As used herein the term "locus group" preferably refers to any combination of two or 
more SNP (single nucleotide polymorphism) loci within approximately 5 cM of each other, 
irrespective of order. 

[0051] As used herein the term "marker" refers to a sequence of DNA that has a specific location 
on a chromosome that can be measured in a laboratory. To be useful, a marker needs to have two 
or more alleles. Common types of markers include, but are not limited to: RELP = restriction 
fragment length polymorphism; SSR = simple sequence repeat (a.k.a. "microsatellite" markers); 
and SNP = single nucleotide polymorphism. Markers may be either located within a known 
gene, or in apparently non-coding regions and not directly associated with a known gene. 
[0052] As used herein the preferred meaning for the term "marker assisted allocation" (MAA) is 
the use of phenotypic and genotypic information to identify animals with superior estimated 
breeding values (EBVs) and the further allocation of those animals to a specific use designed to 
improve the genetic merit of breeding animals for sale or to improve the genetic value of the 
herd. "Allocation" refers to any form of animal management, the selection and distribution of 
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breeding animals, including, but not limited to, marketing the animals as possessing desirable 
characteristics and shipping selected animals to other geographies. In a particular embodiment 
of the invention, "allocation" includes any decision, process, or action that is taken, initiated, or 
considered on the basis of the genetic merit of an animal or animals; where that genetic merit 
was influenced in any way by the genotypic information obtained from the LEPR locus of 
animal(s) or related animals. In this context, the phrase "related animals" refers to the process of 
genotyping an animals and then allocating offspring if/when the offspring's genotype could be 
predicted or assumed, and also covers the use of "allelic peeling" to estimate the genotype of 
ancestors based on the genotype of descendents. 

[0053] As used herein the preferred meaning for the term "Marker assisted selection" (MAS) is 
the use of genotypic information in addition to more traditional phenotypic/pedigree information 
to identify animals with superior estimated breeding values (EBVs) for selection and use as 
breeding animals. 

[0054] As used herein the term "meat color" is used to refer to a the color of uncooked loin and 
ham muscle scored either visually by a trained person using a color scale or objectively using a 
device to measure light reflectance from the cut surface of the meat. One visual industry standard 
is the Japanese Color Score which uses a six point system, with one being the lightest and six the 
darkest. Both Minolta L* and Hunter L* values are often measured objectively using 
instruments manufactured by Minolta Corp. and HunterLab, Inc., respectively. The I* values 
are measures of light reflectance from the cut surface of the meat. Higher L* values correspond 
to higher reflectance and lighter color. 

[0055] As used herein the term "meat quality trait" preferably means any of a group of traits that 
are related to the eating quality (or palatability) of pork. Examples of such traits include, but are 
not limited to muscle pH, purge loss, muscle color, firmness and marbling scores, intramuscular 
fat percentage, and tenderness. 

[0056] As used herein the term "meat tenderness" refers to quantitative evaluation of loin muscle 
tenderness as determined by the Warner-Bratzler shear force test. This measures the force 
required to shear a piece of meat of defined size and orientation that has been cooked to defined 
and controlled specifications. 
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[0057] As used herein the term "muscle mass" refers to the total amount of protein/muscle in the 
animal or carcass. Lean percentage is estimated and the numerator of the formula used is also an 
estimate of the total amount of fat-free lean in the carcass. Estimates of fat-free lean and lean 
percentage are functions of RTUS (real-time ultrasound) measurements of backfat (BF) and loin- 
eye area (LEA), both of which are typically measured in live animals, and body weight. 
[0058] As used herein the term "polymorphism" refers to the variation that exists in the DNA 
sequence for a specific marker or gene. That is, by definition, in order for there to be more than 
one allele for a gene or marker a polymorphism must exist. 

[0059] As used herein a "qualitative trait" is one that has a small number of discreet categories 
of phenotypes. 

[0060] As used herein the term "quantitative trait" is used to denote a trait that is controlled by 
several genes each of small to moderate effect. The observations on quantitative traits are often 
assumed to follow a normal distribution. 

[0061] As used herein the term "quantitative trait locus (QTL)" is used to describe a locus that 
contains polymorphism(s) that has an effect on expression of a quantitative trait. 
[0062] As used herein the term "swine production herd" or "production herd" refers to a 
collection of animals whose primary purpose is to produce pigs that will be shipped to market for 
meat puiposes. 

[0063] As used herein the term "single nucleotide polymorphism (SNP) haplotype" preferably 
refers to a defined combination of SNP alleles from two or more SNP loci on one chromosome. 
[0064] As used herein the term "water holding capacity" preferably refers to measurements of 
drip loss, purge, and cooking loss. The first two measurements are made on uncooked meat and 
measure the ability of the raw meat to hold water. The third measurement refers to cooking loss 
during preparation. It is generally desirable for pork to have low loss of water during storage and 
cooking. 

[0065] According to various aspects of the instant invention the traits to be improved in the pig 
herd may be characterize as either "product quality traits" or "productivity traits". Product 
quality traits for which the instant invention is suitable for improvement include, but are not 
limited to: carcass measurements, meat water holding capacity, meat color, marbling, 
tenderness, and cooking performance. Productivity traits contemplated as part of the instant 
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invention are typically made on live (growing) animals and are known to impact lean growth 
efficiency. Productivity traits contemplated as part of the instant invention include, but are not 
limited to: growth rate, backfat and loin muscle area, feed intake, muscle mass and feed 
efficiency. 

[0066] One of ordinary skill in the art will understand that the product quality traits may be 
measured by a variety of methods. Methods that are known in the art include those that follow. 
Methods for making "carcass measurements" include, but are not limited to: hot carcass weight, 
carcass length, belly thickness, primal weights, Fat-O-Meter fat depth, and Fat-O-Meter loin 
depth. Methods for measuring "water holding capacity" include, but are not limited to: purge 
loss (7 day and 28 day), and drip loss (7 day and 28 day). Methods of evaluating "color" 
include, but are not limited to: determination of Hunter L* (loin), Minolta L* (loin), NPPC 
(National Pork Producers Council) loin color score, and Japanese color score (loin). In addition 
there is a strong correlation between loin and ham pH and measurements of water holding 
capacity and color score. Methods of evaluating "marbling", "tenderness", and cooking 
performance include, but are not limited to: NPPC loin marbling score, NPPC loin firmness 
score, percent intramuscular fat, percent moisture, cooking loss, and Warner-Bratzler shear force 
test. 

[0067] As with "product quality traits", one of ordinary skill in the art will understand that 
"productivity traits" may be evaluated using a variety of techniques which are know in the art. 
Methods for evaluating growth rate include, but are not limited to: evaluating average daily gain 
over various intervals (e.g. birth to 196 days, 90 to 125 days, and 90 to 196 days). Backfat (BF) 
and loin muscle area (LEA) are often measured using real-time ultrasound (RTUS). The 
measurements are typically taken for at least two time points, so as to provide values for BF and 
LEA at specific days of production. Taking measurements at different times also allows for a 
calculation of the change in each of these values over the given interval. Feed intake is often 
measured over a determined interval. For example feed intake may be measured from day 90 to 
day 196. This allows cumulative feed intake over various intervals to be determined (e.g. days 
90 to 196, days 90 to 104, days 90 to 118 and/or etc.) as well as average daily feed intake. Feed 
efficiency is nearly always estimated using measures of growth rate, BF, and LEA, and feed 
intake (if available). 
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[0068] To meet the need for genetic markers associated with desirable porcine traits, the present 
invention provides a method for screening pigs to determine those which will be likely to 
produce offspring with a desirable pork production traits. The method comprises: 1) obtaining a 
sample of genomic DNA from a pig; and 2) analyzing the genomic DNA obtained in 1) to 
determine which leptin receptor allele(s) is/are present and/or determine the animal's status with 
regards to alleles for any markers in linkage disequilibrium with informative markers in LEPR. 
The information collected by this method may then be used in preparing a breeding plan for 
increasing the frequency of the desired allele. 

[0069] The instant invention is further drawn to methods comprising determining which variant 
of the pLEPR polymorphism is extant in an animal, or a plurality of animals, and then using this 
determination to formulate a breeding plan to increase the frequency of the desired allele and/or 
improve the quality of pig offspring produced. 

[0070] It has been discovered that quantitative trait loci (QTLs) for backfat, muscle mass, 
average daily gain (ADG), water holding capacity, meat color, intramuscular fat (marbling), 
tenderness, and cooking loss have been discovered within the area of chromosome 6 where the 
LEPR gene is located (specifically between 100 cM and 130 cM). The inclusion of the pLEPR 
gene (at 122 cM) within the peak of these QTLs is consistent with an association with variation 
in this gene and the variation observed in the described traits {see Example 4). 
[0071] Various embodiments of the current invention provides for the detection of variant alleles 
of one or more single nucleotide polymorphisms within the coding region of the LEPR gene that 
result in, or cause, a polymorphism in one or more amino acid residues of the protein product of 
the pLEPR gene. In preferred aspects of this embodiment of the invention the polymorphisms 
affect pig product quality traits and/or productivity traits. In even more preferred aspects of 
these embodiments the traits are selected from the group comprising, but not limited to feed 
intake, average daily gain, muscle mass, backfat, and water holding capacity, meat color, meat 
pH, intramuscular fat, meat tenderness, and cooking loss. 

[0072] In one aspect of this embodiment of the present invention, the polymorphism comprises 
either a cytosine ("C") or thymine ("T") variant at the nucleotide corresponding to position 609 
of Genbank accession AF1 84172 in the fourth exon of the pLEPR gene. This polymorphism 
produces a pLEPR protein having either a methionine (if the nucleotide is "T") or a threonine (if 



WO 2005/017204 



PCT/US2004/023050 



- 15- 

the nucleotide is "C") at amino acid number 69 of the prepro pLEPR protein. In the animals 
characterized thus far the T variant (containing thymine, encoding methionine) is most 
common {see Example 1). As a shorthand designator, the polymorphism will be referred to as 
"the T69M" polymorphism. 

[0073] Various aspects of these embodiments of the invention provide methods for determining 
the genotype of a particular animal {i.e. genotyping the animal) with respect to the T69M 
polymorphism. That is, determining whether the animal is heterozygous or, alternatively, 
homozygous for one of the variants. These analytical methods include, but are not limited to 
DNA sequencing, primer extension, restriction fragment length polymorphism (RFLP) analysis, 
heteroduplex analysis, single strand conformational polymorphism (SSCP), denaturing gradient 
gel electrophoresis (DGGE), polymerase chain reaction (PCR) both simple and multiplexing (the 
simultaneous amplification of several sequences in a single reaction), real time PCR analysis 
(TAQMAN®), temperature gradient gel electrophoresis (TGGE), allele-specific hybridization, 
oligo-specific hybridization and INVADER® genotyping assays. 

[0074] The current invention provides methods for the testing and/or selection of animals for a 
number of reasons including, but not limited to: breeding (animal husbandry), management, 
forensic purposes, and pedigree analysis. 

[0075] In one embodiment the polymorphism may be identified by an RFLP assay. In one 
aspect of this embodiment the assay may comprise amplifying the pig leptin receptor gene from 
isolated pig genetic material; exposing the gene to a restriction enzyme that yields restriction 
fragments of the gene of varying length. The restriction fragments may then be separated by any 
suitable means. Contemplated methods for the separation of the restriction fragments so as to 
form a restriction pattern, include, but are not limited to such as by gel electrophoresis {e.g., 
using polyacrylamide or agarose gels) or HPLC separation. The resulting restriction fragment 
pattern from the animal is then compared with pig leptin receptor gene that is either known to 
have or not to have the desired marker. If a pig tests positive for the marker, such pig can be 
considered for inclusion in the breeding program. If the pig does not test positive for the marker 
genotype the pig can be culled from the group and used elsewhere. 

[0076] In a preferred embodiment the gene to be analyzed is isolated and replicated using 
oligonucleotide primers and a DNA polymerase to amplify a specific region of the gene that 
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contains the polymorphism. Next the amplified region is digested with a restriction 
endonuclease and the restriction fragments are separated. Visualization of the RFLP pattern is 
by simple staining of the fragments (for example with ethidium bromide), or by labeling either 
the primers or the nucleoside triphosphates used in amplification or both. In a particularly 
preferred aspect of this embodiment the DNA polymerase is a thermostable DNA polymerase 
such as Taq, Pfu, Tfl, or Tli DNA polymerase. 

[0077] Another embodiment of the current invention provides for kits for use in carrying out the 
methods for identifying polymorphisms in the pLEPR gene. These kits comprise the 
components necessary to carry out any method of analysis suitable to detect the polymorphisms 
described herein or known to those of ordinary skill in the art. 

[0078] At a minimum, the kit is a container with one or more reagents that identify a 
polymorphism either in or associated with the pLEPR gene (e.g. in linkage disequilibrium with 
the T69M locus). In one aspect of this embodiment of the invention the kit reagents may 
comprise a set of DNA and/or RNA oligonucleotide primers capable of amplifying a fragment of 
the pLEPR gene that contains the polymorphism. The kit further or alternatively comprise a 
restriction endonuclease enzyme that cleaves the pLEPR gene in at least one place. Other 
possible kit components include, but are not limited to, a DNA polymerase (which may be 
thermostable), a buffer, ribonucleotides and/or deoxyribonucleotides, a reverse transcriptase 
enzyme, and a fluorescent marker. Kits directed to detecting the protein product of the pLEPR 
gene might further comprise a radiomarker and/ one or more antibodies. 

[0079] Other possible components are also considered as part of the instant invention. For 
example, kits comprising components necessary for an immunological assay to detect the allelic 
variant(s) of pLEPR at the protein level. Such components include, but are not limited to, the 
necessary antibodies, buffers, and/or labeling compounds required to perform an enzyme-linked 
immunosorbant assay (BUS A) or any other suitable immunoassay known to those of ordinary 
skill in the art. 

[0080] According to other aspects of this embodiment of the invention the kit may comprise 
oligonucleotide primers suitable for use as a DNA and/or RNA probe or as a primer for DNA 
and/or RNA synthesis. In a preferred aspect of this embodiment of the invention the 
oligonucleotides may comprise or consist of the sequences provided by one or more of the 
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following: SEQ ID NO:l, 2, and 4-9, or the complement of these sequences, or the RNA 
version of these sequences (wherein "U" is substituted for "T"). 

[00S1] Other embodiments of the instant invention provide for oligonucleotides (DNA or RNA) 
that are suitable for use in the kits described above. In preferred aspects of these embodiments 
the oligonucleotides either comprise or consist of the sequences provided in one or more of SEQ 
ID NO:l, 2, and 4-9, or the sequences complementary thereto. 

[0082] Other embodiments of the instant invention provide for methods for producing pigs. 
Generally, the methods comprises analyzing one or more pigs and determining which allele or 
alleles of one or more pLEPR polymorphism(s) each animal possesses. This information 
regarding the allelic composition of the analyzed animal is used as part of a method of managing 
a pig population. 

[0083] In certain aspects of this embodiment of the invention the information collected from the 
analysis of the pLEPR gene in the analyzed pigs in tabulated and utilized, either in isolation or in 
conjunction with other genotypic and/or phenotypic information. In one particular aspect of the 
invention the tabulated information is used as part of a program of marker assisted selection, to 
identify animals with superior estimated breeding values for selection and use as breeding 
animals. 

[0084] In another aspect of this embodiment of the invention the tabulated pLEPR information is 
used as part of a program of marker assisted allocation in order to improve the genetic merit of 
animals to be sold as breeding stock or to improve the genetics of the herd (for example to 
enhance the average estimated breeding value of the herd). 

[0085] In certain aspects of these embodiments of the invention the method includes a breeding 
plan to produce one or more offspring having the desired allelic composition so as to provide the 
qualitative and/or quantitative traits sought. In various aspects of this embodiment the 
information can be used either with or without the assistance of a statistical model/program such 
as BLUP (best linear unbiased prediction) to determine the most effective means to obtain 
animals having the desired traits. In addition to algorithms like BLUP, any other means for 
determining which animals should be bred to each other and/or how the animals should be 
allocated for use in the breeding plan or in the herd, are contemplated by the instant invention. 
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In one particularly preferred aspect of this embodiment of the invention the method is used to 
enhance the accuracy of the estimated breeding value (EBV) for the animals in the herd. 
[0086] Other embodiments of this invention provide for pig herds having more desirable 
characteristics with respect to economic traits selected from, but not limited to, one or more of 
the following: average feed intake and/or average daily weight gain, backfat, muscle mass, 
water holding capacity, meat color, intramuscular fat, meat tenderness, and/or cooking loss. The 
pigs are provided by any of the methods for producing pigs or managing pig populations 
described herein. The instant invention is also drawn to pig offspring produced by any of the 
methods for producing pigs or managing pig populations described herein. 

[0087] Various embodiments of the instant invention are drawn to altering the frequency of a 
pLEPR allele in a selected pig population. The method comprises screening a plurality of pigs to 
identify the nature of an allelic variant in the porcine leptin receptor (pLEPR) gene, wherein said 
allelic variant produces a threonine or methionine polymorphism at amino acid number 69 of the 
prepro-pLEPR protein. Such screening can be accomplished either by directly determining the 
DNA sequence or by any other suitable method, for example by determining the sequence of the 
pUEPR gene product protein or by identifying either a SNP or SNP haplotype known to be in 
linkage disequilibrium with a particular allelic variant. Once the nature of the polymorphism is 
known for each animal, then those pigs having the desired allelic makeup are selected. These 
pigs are then allocated for use according to a breeding plan designed to achieve the desired 
change in the pig populations allelic frequency. This breeding plan may be designed to increase 
the frequency of a particular allelic profile. One aspect of this embodiment of the invention 
includes employment of a plan to fix the allele in a given pig population. Alternatively, the 
breeding or managerial plan may be designed to decrease the allele or to provide for a more 
"balanced" occurrence of the allele. 

[0088] Another embodiment of the instant invention provides a method for enhancing meat 
production (that is improving either the quality or the quantity of the meat) in a herd by a method 
comprising modifying the herd genetics (e.g. the frequency of a particular pLEPR gene allele) as 
provided herein. One aspect of this embodiment of the invention provides for a method wherein 
the EBV of the herd is improved over time with respect to the trait of meat production by, for 
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example, increasing the frequency of the LEPR gene allele which has been shown to be linked to 
increased meat production. 

[0089] Another embodiment of the invention provides for pig offspring produced using the 
methods and/or kits described herein. In various aspects of this embodiment individual offspring 
or litters; multiple offspring or litters; or entire herds may be produced using the methods 
described. That is, these offspring may be produced by methods comprising identifying animals 
having a desired pLEPR gene polymorphism and using these animals in a breeding program to 
produce offspring. 

[0090] The instant inventors have found an association between the T69M locus and various 
growth-related phenotypes (see Example 4). An analysis of 2625 pigs from a single commercial 
line, showed that the presence of the "C" allele had a statistically significant correlation with a 
positive effect on: early ADG (average daily gain from day 0 to day 90 of life); late ADG 
(average daily gain from day 90 to day 165 of life), loin muscle pH, and loin muscle color, and 
drip loss. There was a small negative effect of the "C" allele on backfat, Le. backfat was slightly 
increased. 

[0091] In addition, ninety seven (97) SNP markers, representing 38 loci on porcine chromosome 
6 (SSC6) were genotyped on a panel of 1,442 pure line pigs from the same commercial line. The 
loci selected for SNP discovery were spread across an approximately SO cM region on SSC6 
which included the LEPR locus and the SNP producing the T69M mutation. Linkage 
disequilibrium analysis was used to identify both individual SNPs and SNP locus groups (for up 
to three adjacent/nonadjacent SNPs that locate within 1 centiMorgan (cM)) that were 
significantly associated with growth-related phenotypes (i.e. backfat thickness, leanness, off-test 
weight and weight gain). All 97 SNPs and possible locus combinations of two and three SNP 
located within 1 cM were assessed for association with all phenotypes. At least four SNPs (plus 
several locus groups containing these SNPs) were found to be significantly associated with 
backfat thickness, corrected for either age or weight. One of these SNPs included T69M and the 
other three mapped within 3 cM of T69M as estimated by linkage analysis (see Table 5). 
[0092] One significant feature of the present invention is that it is drawn to a pLEPR gene 
polymorphism in the coding region of the gene and causes an amino acid change in the protein 
product of the pLEPR gene. Li contrast, previously published polymorphisms are not 
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characterized as causing an amino acid change, instead they are believed to occur in the intronic, 
or non-coding, portion of the gene. 

[0093] Other aspects of the instant invention provide for methods of identifying one or more 
single nucleotide polymorphism(s) in linkage disequilibrium with the T69M polymorphism as 
described in the EXAMPLES, below. Briefly, various aspects of this embodiment of the 
invention comprise identifying at least one large-insert genomic clone containing all or a portion 
of the pLEPR gene. This large-insert genomic clone may be obtained from a porcine genomic 
library and may be in any suitable format. Suitable formats include, but are not limited to, a 
bacterial artificial chromosome (BAC), a yeast artificial chromosome (YAC), PI, a cosmid, a 
fosmid, a phage, and a plasmid. 

[0094] Next large-insert clones containing all or part of the pLEPR gene are identified by any 
suitable means. According to one aspect of this embodiment the clones are identified by 
hybridization with a DNA or RNA probe comprising all or a portion of the pLEPR gene. Once 
one or more clones comprising all or a portion of the pLEPR gene are identified then the 
sequence of all or a portion of such clones may be determined. In addition, sequences 
representing genes and expressed sequence tags (EST) and markers (e.g. microsatellites) that 
have been placed on either physical (e.g. radiation hybrid) or linkage maps and that appear to be 
in proximity to the pLEPR gene are identified and used to select BAC clones containing these 
sequences. Furthermore, sequences selected from comparative maps (e.g. human-porcine) that 
appear to be in the proximity of pLEPR can also be used to screen and select BAC clones. 
Following the sequencing of these clones, portions of the clones comprising regions in close 
proximity to the pLEPR gene can be identified, referred to herein as "target regions." In this 
context, "close proximity" refers to any chromosomal distance over which linkage 
disequilibrium may exist, preferably up to 5 cM (roughly equivalent to 5 million base pairs). 
Factors influencing linkage disequilibrium vary between populations and include effective 
population size, mating structure, generation interval, ancestry, and other factors. Once one or 
more target regions are identified a panel of animals is screened to determine the sequence of 
their genomes in the areas corresponding to the target regions. The data generated from this 
screening is then analyzed to identify any single nucleotide polymorphisms (SNPs) present 
therein. The nature of the T69M allelic variant is also determined for each of these animals. 
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FinaUy, the SNP data is analyzed with respect to each newly identified SNP to determine which 
of the newly identified SNPs is in linkage disequilibrium with the T69M polymorphism. 
[0095] According to various aspects of this embodiment of the invention, SNPs identified as 
being in linkage disequilibrium with the T69M polymorphism are useful as markers for use in 
any of the methods described herein. 

DESCRIPTION OF THE SEQUENCE LISTINGS 
[0096] The following sequence listings form part of the present specification and are included to 
further demonstrate certain aspects of the present invention. The invention may be better 
understood by reference to one or more of these sequences in combination with the detailed 



description of specific embodiments presented herein. 





j Descnption 


1 

1 


LiirK-KrLr-r i primer 


9 


l^i^r iv-xCrLvJr -K l primer 


D 


urenoank accession no: AJr 1 o4173 


4 


LEPR-311-F 


5 


LEPR731-R 


6 


Forward primer for T69M TAQMAN® assay 


7 


Reverse pnmer for T69M TAQMAN® assay 


8 


Probe for T69M TAQMAN® assay 


9 


Probe for T69M TAQMAN® assay — 


10 


Sequence from LEPR exon 


11 


LEPR sequence used to select BAC 335C21 


12 


Sequence from BAC 335C21, which was used to identify informative SNPs 


13 


Sequence used to select BAC 036M15, which in close proximity to LEPR locus 


14 


Sequence from BAC 036M15, which was used to identify informative SNPs 


15 


Sequence used to select BAC 069P03 in close proximity to LEPR locus 


16 


Sequence from BAC 069P03, which was used to identify informative SNPs 


17 


Forward primer for TAQMAN® assay number 183482 


18 


Reverse pnmer for TAQMAN® assay number 183482 


19 


vicProbe for TAQMAN® assay number 183482 


20 


famProbe for TAQMAN® assay number 183482 


21 


Forward pnmer for TAQMAN® assay number 180851 


22 


Reverse pnmer for TAQMAN® assay number 1 8085 1 i 


23 


vicProbe for TAQMAN® assay number 180851 


24 


famProbe for TAQMAN® assay number 1S0S51 


25 


Forward pnmer for TAQMAN® assay number 182553 


26 


Reverse primer for TAQMAN® assay number 1 82553 


27 


vicProbe for TAQMAN® assay number 1 82553 j 


28 


famProbe for TAQMAN® assay number 1 S2553 


29 


Forward genomic pnmer derived from BAC clone 069P03 


30 


Reverse genomic primer derived from BAC clone 069P03 H 
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31 


Forward genomic primer derived from BAC clone 036M15 


32 


Reverse genomic primer derived from BAC clone 036M15 


33 


Forward genomic primer derived from BAC clone 335C21 


34 


Reverse genomic primer derived from BAC clone 335C21 


35 


Forward primer to screen BAC library (705625F) 


36 


Reverse primer to screen BAC library(705625R) 


37 


Forward primer to screen BAC library (AR024A1 IF) 


38 


Reverse primer to screen BAC library (AR024A 1 1R) 


39 


Forward primer to screen BAC library (3661588F) 


40 


Reverse primer to screen BAC library(3661588R) 


41 


Amplicon amplified using primers derived from BAC clone 069P03 


42 


Amplicon amplified using primers derived from BAC clone 036M15 


43 


Amplicon amplified using primers derived from BAC clone 335C21 i 



[0097] The following examples are included to demonstrate preferred embodiments of the 
invention. It will be appreciated by those of skill in the art that the techniques disclosed in the 
examples which follow represent techniques determined by the inventors to function well in the 
practice of the invention, and thus can be considered to constitute preferred modes for its 
practice. However, those of skill in the art should, in light of the present disclosure, appreciate 
that many changes can be made in the specific embodiments which are disclosed and still obtain 
a like or similar result without departing from the scope of the invention. 

EXAMPLES 

Example 1: Characterization of the "C"/"T" polymorphism in swine 

[0098] PCR primers, corresponding to SEQ ED NO:4 and SEQ ID NO:5, were designed to 
amplify a portion of the LEPR gene containing a small intron and the beginning of the coding 
region. These primers were used as part of a PCR reaction using the DNA from IS animals as 
template. The resultant amplicons were sequenced and analyzed for polymorphisms. This 
analysis resulted in the identification of a polymorphism at nucleotide 299 of the amplicon 
(corresponding to position 609 of Genbank accession API 84173, SEQ ID NO:3). 
[0099] The set of 18 animals used for polymorphism discovery were from a first ("Line A"; 
Pietrain) and a second ("Line B"; Duroc) commercial line. DNA was extracted from either ear 
or tail tissue using commercially available DNA extraction materials (Qiagen N.V., Venlo, 
Netherlands). DNA was subjected to PCR amplification using oligonucleotide primers SEQ ID 
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NO: 4 and SEQ ID NO: 5. Amplified fragments were sequenced in both directions using the 
amplification primers. Resultant DNA sequences were called, aligned, and characterized for 
polymorphism using the Phred/Phrap/Consed/Polyphred software package developed and 
distributed by Phil Green (University of Washington, Seattle, Washington). The polymorphism 
described herein was detected at position 299 of the amplified sequence, and was determined to 
alter the amino acid sequence of the LEPR protein. 
Example 2: Discovery of SNPs in close proximity to LEPR 

[0100] PCR primers were designed from a portion of the complete coding sequence of pLEPR 
(SEQ ID NO: 11) from a porcine EST sequence that was in close proximity to pLEPR on a 
porcine radiation hybrid map (SEQ ID NO: 13) and from a porcine EST sequence that was 
homologous to a sequence obtained from the human sequence map and in close proximity to 
human LEPR (SEQ ID NO: 15). These PCR primers (SEQ ID NO:35-^0) were used to screen 
BAC clones from a porcine bacterial artificial chromosome library (RPCI-44 see 
bacpac.chori.org/mporcine44.htm) and select a clone representing each sequence/locus (BAC 
clones 069P03, 036M15, and 335C21, respectively). These three BAC clones were then 
subcloned and approximately 48 subclones were randomly selected and sequenced. High quality 
subclone sequences that did not contain known porcine repetitive elements were then selected for 
another round of primer design. Due to the fact that only partial and often times non overlapping 
sequence was obtained for each of the selected BAC clones, the sequences selected for primer 
design usually did not include any of the original sequence (SEQ ID NO: 11, 13, and 15) used to 
screen the BAC library. 

[0101] The genomic sequences derived from BAC clones 335C21, 036M15 and 069P03 and 
were used to identify SNPs in LD with the T69M polymorphism are represented by SEQ ID 
NOs: 12, 14, and 16, respectively. The PCR primers that were designed from these sequences 
and that were then used to amplify genomic DNA template from a panel of 18 animals are shown 
as SEQ ID NO:29-34. The sequences of the amplicons produced using these primers, together 
with the location of polymorphisms identified by aligning and comparing the sequences from 
each of the 18 animals used in the discovery panel are provided as SEQ ID NO:41-^3. 
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[0102] The TAQMAN® SNP assay developed for the SNP identified in SEQ ID Nos:41, 42, and 
43, were assigned assay numbers 183482, 180851, and 182553, respectively. The details for the 
primers and probes used these TAQMAN® assays are provided in tables 1-3. 



fwd Primer 


5'-GGCAGCTGTAACTGGTTACGAA-3' (SEQ ID NO: 17) 


rev Primer 


5'-TCGC AGCTC ATATTGAATAACGATGT-3 ' (SEQ ID NO:18) 


vie Probe 


5'- AAGTTCCA AATACTCTTTC-3 ' (SEQ ID NO: 19) 


vicAllele 


G 


famProbe 


5-AAGTTCCAAATACTATTTC-3 ' (SEQ ID NO:20) 


famAllele 


T 

1 t 



Ta ble 2: Probes and primers for TAQMAN© assay number 180851 



fwd Primer 


5'-CAGACCCTCTGATATTTGGAAAAGCA-3 * (SEQ ID NO:21) 


rev Primer 


5-GCC AGG AT AATC ATTTGAGT AT AAG A A A AG AAC-3 ' (SEQ 
ID NO:22) 


vicProbe 


S'-ACAGGAGCTACTAAAAT-S' (SEQ ID NO:23) 


vicAllele 


C 


famProbe 


5-C AGGAGCT ATTAAAAT-3 ' (SEQIDNO:24) 


famAllele 


T 



Ta ble 3: Probes and primers for TAQMAN® assay number 182533 



fwd Primer 


5-ACATTCTAAGACAACCGAAATGGCA-3' (SEQ ID NO:25) 


rev Primer 


5'-CT AGGGATCT ATTTTTC ACTTTTGT AAGTTC ATT-3 ' (SEQ ID 
NO:26) 


vicProbe 


5'- ATAATTTTCATAAAGACCCACTAAT-3 * (SEQ ID NO:27) 


vicAllele 


A 


famProbe 


5'-C ATAAAGGCCC ACTA AT-3 * (SEQK>NO:28) 


famAUele 


G 



Example 3; Methods for Genotyping the T69M locus 

[0103] Several methods exist to determine allelic composition at the LEPR T69M 
polymorphism. Such methods include, but are not limited to, PCR amplification and sequencing 
using SEQ ID NO: 4 and SEQ ID NO: 5 or other suitable primer pairs consisting of DNA 
sequence flanking the polymorphism, RFLP analysis using amplification primers SEQ ID NO: 1 
and SEQ ID NO: 2 or other suitable primers flanking the polymorphism in conjunction with a 
restriction endonuclease such as BsrDl or other suitable enzyme to discriminate between the "C" 
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and "T" alleles in the amplified DNA, real time PCR analyses (TAQMAN ®) involving DNA 
amplification and probe hybridization where the hybridization probes are labeled and 
discriminate between the allelic forms, and other methods readily performed by those skilled in 
the art (see Table 4). 



Table 4 



Fwd Primer 


5*-TTCAACTTTGAATGGACATGATGAG-3 > (SEQ ID NO:6) 


rev Primer 


5 , -GTGGAAAGTTGTTTTAGAAGATAAGTTTGA-3 , (SEQ ID NO:7) 


vicProbe 


5 '-TGTTGAAACGGAACTT-3 * (SEQ ID NO:8) 


vicAllele 


C 


famProbe 


5 * -TGTTG A AATGG AACTT A-3 ' (SEQ ID NO:9) 


famAllele 


T 



Example 4: Associat ion between the pLEPR (T69M) polymorphism with production- and 
meat quality-related traits 

[0104] To carry out the association analysis between the pLEPR (T69M) polymorphisms and 
production and meat quality traits, almost 3,000 "Line A" pigs, representing more than 100 
paternal half-sib families, with a combination of production and meat quality trait records were 
genotyped. 

[0105] Analysis was performed via the following steps. First, biologically impossible data and 
the phenotypes that were outside the range of mean plus/minus 4 phenotypic standard deviations 
were excluded. The phenotypes were then modeled to account for the effects of season-year- 
farm-building, sex, and age effect, and were excluded if their residuals were outside the range of 
residual mean plus/minus 4 residual standard deviations. Remaining phenotypes for all animals 
were pre-adjusted using the following two models: 

Phenotype = cgp + sex + age +sire + residual [1] 
Phenotype = cgp + sex + age +sire + dam + residual [2] 

Where phenotype in model [1] denotes phenotype for daily gain of body weight for the whole 
growth period (WD A), daily gain of body weight for the test period (ADG2), backfat depth (BF) 
and loin eye area (LEA) measured at the 10 th rib on the day of off test, and phenotype in model 
[2] denotes daily gain of body weight for the period from birth to on test (ADG1); cgp represents 
contemporary group that was formed to account for season-year-farm-building effect; sire and 
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dam are class variable to account for sire and dam effect. The association between the candidate 

gene polymorphism and pre-adjusted trait phenotypes vas evaluated using following model: 

Phenotype = genotype + residual 
Phenotype = # allele + residual 

The effect of genotype or allele were estimated via a least squares procedure, and the p value was 
estimated based on an F distribution with the degrees of freedom equal to number of genotypes - 
1 or number of alleles -1, respectively. 

[0106] Not all animals had a complete set of phenotypic records, thus in Tables 5 and 6 the 
number of animals included in the analysis is indicated. 

Table 5. Deviation from adjusted population mean for the three genotypic classes of 
pLEPR T69M for production traits 





GENOTYPE 






Trait 


C/C 


C/T 


T/T 


p-value 


# animals 


On test weight (lbs) 


7.09 


3.51 


1.72 


0.0039 


1056 


Off test backfat 
(inches) 


0.0109 


0.0084 


-0.0046 


0.0012 


2590 


Off test lean percent 


-0.272 


-0.149 


0.059 


0.0242 


2589 



Table 6. Deviation from adjusted population mean for the three genotypic classes of 
pLEPR T69M for meat quality traits 





GENOTYPE 






Trait 


C/C 


C/T 


T/T 


p-value 


# animals 


Loin Color 
(Hunter L* @ day 7) 


-0.580 


-0.424 


0.008 


0.0158 


388 


Loin pH @ day 7 


0.032 


0.011 


-0.013 


0.0315 


460 



[0107] The predicted impact on average phenotypic values for weight gain, backfat, loin color, 
and loin pH by fixing the "C" allele in Line A boars was estimated {see Table 7). The initial 
frequency of the "C" allele in this boar line is 22% and changes in pure line offspring phenotypic 
values were estimated assuming frequency of the "C" allele were increased to 100%. 
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Table 7. Impact of fixing "C" allele in terminal boar line on pure line offspring 



Trait 


i Predicted change 


Percent of mean 


P-value 


Gain on test 
(lbs from 0-90 days) 
(lbs from 90-165 days) 


+1.48 
+2.01 


1.4 
1.4 


<0.05 
<0.001 


Backfat 
(inches @ day 165) 


+0.01 


2.3 


<0.001 


Loin color 
(Hunter L* @ day 7) 


-0.375 


1.0 


<0.05 


Loin pH @ day 7 


+0.021 


0.36 


<0.01 



[0108] Ninety seven (97) SNP markers, representing 38 loci on porcine chromosome 6 (SSC6) 
that were genotyped on a panel of 1,442 pure line pigs. Analysis was performed via the 
following steps: first, biologically impossible data and the phenotypes that were outside the 
range of mean plus/minus 4 phenotypic standard deviations were excluded. Next, phenotypes 
were modeled to account for the effects of season-year-farm-building, sex, and age effect, and 
were excluded if their residuals were outside the range of residual mean plus/minus 4 residual 
standard deviations. Remaining phenotypes for all animals of the entire pure line were pre- 
adjusted using the following two models: 
Phenotype = cgp + sex + offage + residual 
Phenotype = cgp + sex + offwt + residual 

where cgp was designed to account for season-herd-building effect, offage and offwt denote age 
and body weight at off test (i.e., when measurements are taken), respectively. 
[0109] The second step was to form locus groups, and to calculate probabilities for each possible 
haplotype pair for each animal. As a preparation, a linkage map for these 38 loci was 
constructed using genotype information of approximately 3000 animals, in combination with 
their radioactive hybridization information. Within each locus, the order of SNP was then 
arbitrarily assigned, and the linkage distance between adjacent within locus SNP was assumed to 
be 0.01 centiMorgan. Based on their linkage map information, every combination of 1 to 3 SNP 
markers that locate within the distance of 1 centiMorgan forms a locus group. For all 97 linked 
SNP, the probability of likely linkage phase of sires with SNP genotypes were calculated 
conditional on pedigree and the SNP genotypes of their parental, mate, and progeny genotype 
information, using a very efficient algorithm. Conditional on sire linkage phases, probability of 
each possible haplotype pair was calculated for each animal and for each locus group. 
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[0110] The third step was to evaluate the association between preadjusted trait phenotypes and 
haplotype pairs of animals using following model: 

K 

Yi = ZPkXjk+e; 

k=l 

where yj and ej are the preadjusted trait phenotype and the residue for animal i, respectively; \ lk 
denote the sum of probability of both the paternal and maternal haplotype being k, 0 k is 
regression coefficient for haplotype k, and K is the total number of haplotypes in the population. 
The Type I error rate ( p value) was estimated by performing 50,000 random permutations of 
phenotypes among paternal half-sibs. The results showed that at least four SNP markers (very 
tightly linked to the LEPR locus) and four locus groups that were significantly associated with 
backfat thickness. Table 8 shows the F-statistic, p value, linkage map position, frequency of the 
favorable allele and the estimated effect of fixing the favorable allele for these four SNPs. Table 
9 shows the equivalent information for the four SNP haplotypes found to be significantly 
associated with backfat thickness. In the case of each of these SNP combinations, two haplotypes 
accounted for >99% of the observed genotypes, thus these haplotypes were essentially biallelic. 



Table 8. SNPs significantly associated with backfat thickness on SSC6 



SNP/Assay # 


Map position 
(cM) 


F-statistic 


p-value 


Favorable 

allele 

frequency 


Fixation 
effect 


180851 


132.3 


20.57 


<0.0001 


0.719 


-0.00438 


182553 


133.3 


21.18 


<0.0001 


0.720 


-0.00443 


LEPR(T69M) 


133.3 


20.67 j 


<0.0001 


0.721 


-0.00436 


183482 


135.9 


13.85 


0.0001 


0.486 


-0.0057 



Table 9. SNP haplotypes significantly associated with backfat thickness on SSC6 



Locus 
Group ID 
No. 


SNPs comprising 
haplotype 


F-statistic 


p-value 


Favorable 
haplotype 
frequency 


Fixation 
effect 


125 


180851 + 182553 


21.55 


<0.0001 


0.870 


-0.00204 


126 


180851 + 182553 
+ LEPR(T69M) 


21.28 


<0.0001 


0.872 


-0.00200 


127 


180851 + 
LEPR(T69M) 


21.34 


<0.0001 


0.872 


-0.00200 


129 


182553 + 
LEPR(T69M) 


21.28 


<0.0001 


0.873 


-0.00201 
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Additional evidence to support the contention that LEPR polymorphisms are associated 
with productivity and meat quality traits. 

[0111] A resource population for QTL discovery was created by crossing Pietrain boars with 
Duroc sows. The Fl generation was intercrossed to produce a F2 generation in which alleles 
differing between the two founder lines were expected to be segregating. In total, 1,600 F2 
progeny were generated and specific productivity and meat quality phenotypes were measured 
on at least 1,000 of these animals (depending on the specific trait). Approximately half of the F2 
generation and their parents and grandparents were genotyped for 135 microsatellite markers 
spaced across all 18 autosomes, autosomes (9 microsatellite markers from SSC6). 
[0112] Described below is a typical analysis procedure. First, performance trait phenotypes for 
all F2 animals were pre-adjusted using the following two models: 
Growth phenotype = cgp + sex + age + residual 

RTUS phenotype = cgp + sex + offwt + residual 

where cgp was designed to account for season-herd-building-pen effect, and formation of cgp 
was different for different ages; age denotes the age when the measurement was taken; RTUS 
denotes real time ultrasound measurements of backfat or loin-eye area. For meat quality traits, 
phenotypes were pre-adjusted using following models: 
MQ1 = sdate + residual 

MQ2 = sdate + sex + residual 

MQ3 = sdate + sex + sage + residual 

where MQ1 denotes drip loss after 7 days, drip loss after 28 days, or purge loss after 28 days; 
MQ2 denotes intramuscular fat, marbling score, or percent moisture; MQ3 denotes Warner- 
Bratzler shear force; sdate is slaughter date fitted as a class variable, and sage is slaughter age 
fitted as a covariate. 

[0113] Second, grandparental line origin of F2 offspring was traced for each marker, as 
described by Haley et al. (1994). Two coefficients (x a and x d ) were calculated as the difference 
in probability of being homozygous and the probability of being heterozygous. 
[0114] Third, the association between the preadjusted phenotype was evaluated as: 
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where yi is the preadjusted trait phenotype for animal i, 0 O , 0 a , and (3 d are regression coefficients. 
For each putative QTL position, residual sum squares is minimized to estimate regression 
coefficients and test statistic, and a chromosome was searched in incremental 1 cM steps, by 
performing analysis for each putative QTL position. Chromosomewise Type I error rate (p value) 
was estimated by performing 10,000 to 20,000+ random permutations of F2 phenotypes within 
each paternal half-sib family to determine empirically the proportion of times observed test 
statistics occurred by chance. 

[0115] Porcine LEPR is closely linked (approximately 1 cM on MARC map) to SW1881, which 
was located at 121 cM on the linkage map constructed for the markers genotyped on SSC6. 
Tables 10 and 11 list traits that had a' significant (p <0.05) F-statistic at 122 cM for productivity 
and meat quality traits, respectively. 

Table 10. F-statistic and probability for productivity traits that had QTL coinciding with 
the predicted location of pLEPR 



TRAIT 


F-STATISTIC 


/7-Value 


Average daily gain from d 0 to 56 


13.21 


<0.0001 


Average daily gain from d 21 to 56 


10.88 


<0.0001 


Average daily feed intake from d 90 to 196 


12.31 


<0.0001 


Backfat at d 90 


60.72 


<0.0001 


Backfat at d 124 


73.40 


<0.0001 


Backfat at d 160 


81.34 


<0.0001 


Backfat at d 196 


75.90 


<0.0001 


Change in backfat from d 90 to 196 


48.13 


<0.0001 


Change in loin eye area from d 125 to 196 


7.73 


0.0052 
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Table 11. F-statistic and probability for meat quality traits that had QTL coinciding with 
the predicted location of pLEPR 



TRAIT 


F-STATISTIC 


p-Value 


Drip loss after 7 days 


6.79 


0.0125 


Drip loss after 28 days 


5.59 


0.0347 


Purge loss after 28 days 


6.77 


0.0122 


Intramuscular fat 


31.45 


<0.0001 


Marbling score 


22.69 


<0.0001 


Percent moisture 


8.23 


0.0047 


Wamer-Bratzler shear force 


12.53 


<0.0001 



[0116] Discovery of QTL for average daily gain, backfat and loin eye area in the F2 resource 
population closely associated with the predicted location of pLEPR on chromosome 6 
corroborates the association between polymorphisms in pLEPR and growth rate, backfat and 
leanness discovered using "Line A" pigs. In addition, discovery of QTL for drip loss and purge 
loss in the F2 resource population closely associated with the predicted location of pLEPR on 
chromosome 6 also supports the association between pLEPR polymorphisms and loin color and 
pH. Although a QTL of large effect for loin pH and color score on chromosome 6 was not 
specifically identified by these particular analyses, these two traits are known to be highly 
correlated with drip loss and purge. 

[01171 Genetic correlations between drip loss and pH, drip loss and color, and pH and color for 
data collected from our F2 animals were -0.66, 0.45, and -0.37, respectively. In addition, Huff- 
Lonergan et dL % 2002, reported significant (p <0.0001) phenotypic correlations between percent 
drip loss and pH and loin color score. These correlations are consistent with the observation that 
meat with a higher pH tends to be darker in color (lower L* values) and have less drip loss. 
Thus, in view of the information provided herein, it would be expected by one of skill in the art 
that polymorphisms in pLEPR would also be associated with drip loss and purge as well as loin 
color score and pH. 

[0118] The association between pLEPR polymorphisms and intramuscular fat and meat 
tenderness did not exceed statistical significance, primarily because these phenotypes were not 
measured for the majority of animals from which DNA was taken for pLEPR T69M genotyping. 
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However, backfat and percent intramuscular fat are highly positively correlated traits. Huff- 
Lonergan et al y 2002, reported a significant association (0.45, p <0.0001) between 10 th rib 
backfat and percent lipid. In addition, Ovilo et al y 2000 and de Koning et al., 1999 also detected 
QTL for intramuscular fat on porcine chromosome 6. Therefore, based on the results provided 
herein, one of ordinary skill in the art would not find it surprising that the work disclosed herein 
also identifies QTL for intramuscular fat in the same location. 

[0119] Genetic correlations between intramuscular fat and percent moisture and intramuscular 
fat and Warner-Bratzler shear force for data collected from our F2 animals were -0.83 and -0.63, 
respectively. These correlations are consistent with the observation that increased intramuscular 
fat decreases muscle protein and associated water and that meat tenderness increases (shear force 
decreases) with increased marbling or intramuscular fat. Thus, it is reasonable to assert that 
polymorphisms in pLEPR are associated with measurements of intramuscular fat, moisture and 
tenderness. 

Example 5: Identification of Single Nucleotide Polymorphism^) (SNP(s)) in linkage 
disequilibrium with the pLEPR T69M SNP 

[0120] A person skilled in the art could discover Single Nucleotide Polymorphisms (SNPs) in 
Linkage Disequilibrium (LD) with pLEPR T69M by processes similar to (but not limited to) the 
following: 

[0121] The skilled artisan could identify a large-insert clone containing either the pLEPR gene 
or sequences in close proximity. Such a clone could be a bacterial artificial chromosome (BAC), 
yeast artificial chromosome (YAC), PI phage, cosmid, fosmid, phage, or plasmid constructs. 
Obtaining this clone could involve hybridization to genetic libraries with labeled DNA or RNA 
probe, or by iterative PCR, using primers and/or probes known to amplify sequences at or near 
the pLEPR gene. 

[0122] The molecule containing the pLEPR gene could then either be sequenced directly or be 
subcloned and then sequenced to identify specific DNA sequences known to exist in close 
proximity to (or flanking) the pLEPR gene. For the purposes of this example, these flanking 
sequences are referred to as "target" sequences. The number of target sequences obtained is 
relevant insofar as the presence of more target sequences proportionally increases the likelihood 
of identifying a SNP in LD with T69M. 
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[0123] Once target sequence is identified, primers suitable for use in Polymerase Chain Reaction 
(PCR) amplification of target DNA from a panel of animals (the "SNP discovery panel") could 
be designed. Target DNA derived from the SNP discovery panel could then be sequenced and 
any SNPs present in the discovery panel that are in LD with the T69M could be identified, if 
present. By definition, these SNPs would be physically located in close proximity to the pLEPR 
gene. The number of SNPs thus identified is relevant insofar as the discovery of more SNPs 
proportionally increases the likelihood of identifying a SNP in LD with T69M. 
[0124] Once a set of SNPs to be tested for LD has been obtained, the skilled artisan could 
conduct experiments to calculate LD between T69M and the individual SNPs identified in the 
set. These experiments are conducted by identifying an independent panel of animals (the "LD 
panel") that are unrelated and representative of as many phylogenetically distinct breeds as 
practicable. 

[0125] All animals within the LD panel could be genotyped for all SNPs within the set, as well 
as for the T69M SNP. Each genotype represents two alleles, one each from a distinct 
chromosome, one maternal and one paternal. Therefore, a "phase" relationship can exist 
between alleles at two loci. For example, assume two genes (A and B), each have two alleles 
(Al and A2; Bl and B2), thus there are four combinations of alleles for each chromosome (Al- 
Bl, A1-B2, A2-B1, and A2-B2). Most genotypes can be deconstructed to derive the two 
component chromosomes. 

[0126] Linkage disequilibrium (LD) is then measured between two loci by using allele frequency 
data from the LD panel to calculate the expected frequency of the four combinations of alleles 
for each chromosome (Al-Bl, A1-B2, A2-B1, and A2-B2) assuming random distribution 
(frequency of Al-Bl = frequency of Al x frequency of Bl, etc.). These frequencies are 
compared to the observed distribution of allele combinations. If the frequencies are significantly 
different (therefore non-random), the two loci are said to be in LD. 

[0127] If two alternate combinations of alleles are never observed (for example, A1-B2 and A2- 
Bl), then the two loci are said to be in complete LD . This situation provides the perfect 
opportunity to use one locus as a proxy for the second locus when genotyping animals. 
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